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Abstract  —  The  scattering  coefficient  of  seawater  as  a 
function  of  concentration  of  hydrosol  particles  is  calculated. 
The  approach  used  is  based  on  the  Maxwell's  equations  in  a 
stochastically  scattering  seawater.  The  water  is  modeled  as 
thermally  fluctuated  medium  filled  with  the  hydrosol 
particles.  It  is  found  that  the  scattering  coefficient 
quadratically  depends  on  concentration  when  the  concentration 
of  scatterers  is  very  small.  The  scattering  coefficient  is  linear 
to  concentration  of  scatterers  at  values  typical  to  the  open 
ocean.  At  the  values  of  concentrations  typical  to  coastal 
waters  the  dependence  on  concentration  weakens  and  reaches 
the  saturation  at  very  high  values. 

INTRODUCTION 

In  a  majority  of  publications  available  today  the  light 
scattering  coefficient  by  seawater  is  considered  as  linearly 
dependent  on  concentrations  of  hydrosols.  The  experiments  by 
Prieur  and  Sathyendranath  [1]  show  that  at  certain 
concentrations  of  chlorophyll  C^-,  typical  to  coastal  waters, 
the  dependence  of  seawater  absorption  coefficient  on  C<-  is 
nonlinear.  Clark  and  Backer  [2]  showed  that  the  scattering 
coefficient  of  seawater  is  non-linearly  dependent  on 
chlorophyll  concentration  that  is  strongly  correlated  to  the 
concentration  of  scattering  matter  of  biologic  origin. 

In  this  paper  an  attempt  is  made  to  develop  an  approach 
to  calculate  the  scattering  coefficient  of  seawater  as  a  function 
of  concentration  of  hydrosol  particles  C^.  The  approach  is 
based  on  the  solutions  of  the  Maxwell's  equations  in  a 
stochastically  scattering  medium  (seawater).  The  seawater  is 
modeled  as  thermally  fluctuated  medium  filled  with  the 
hydrosol  scatterers. 

The  final  result  of  this  paper  is  Eqn.  (22)  for  the  seawater 
scattering  coefficient.  The  scattering  coefficient  b  is  linear 
with  concentration  of  scattering  particles  at  values  typical  to 
the  open  ocean.  The  coefficient  quadratically  depends  on 
concentration  when  the  concentration  is  very  small  (typical 
for  the  Sargasso  Sea  waters).  At  the  concentrations  close  to 
the  values  that  are  typical  to  coastal  waters,  the  dependence  on 
concentrations  weakens  and  reaches  saturation  at  very  high 
values.  The  results  of  this  paper  can  explain  some 
experimental  data  obtained  in  turbid  coastal  waters.  They  also 
can  give  a  reasonable  explanation  to  the  phenomenon  that 
some  very  clean  ocean  waters  seem  to  be  more  transparent 
than  the  distilled  water. 

APPROACH 


mathematical  formalism  of  scattering  in  stochastic  medium  is 
identical  to  the  formalism  of  quantum  statistical  mechanics 
[3].  The  photons  themselves  are  always  quantum  particles. 
For  these  reasons  we  have  chosen  to  use  the  quantum- 
mechanical  statistical  approach  formulated  in  Ref.  [4]. 

In  this  paper  the  scattering  coefficient  on  hydrosol 
particles  is  calculated  through  the  dielectric  permittivity  of  the 
hydrosol  component  of  the  water.  The  dielectric  permittivity 
(DP)  is  a  constituent  part  of  the  Fourier  transform  of  the 
Green's  function  of  the  Maxwell's  equations  in  seawater.  By 
definition,  the  Green's  function  is  a  solution  of  these 
equations  when  a  source  function  is  assumed  as  an  infinitely 
short  and  localized  at  one  point  light  pulse  [5].  The  approach 
adopted  here  is  based  on  the  theory  proposed  by  the  author  in 
Refs.  [6,  7]  but  abandoned  at  the  time  due  to  the  lack  of 
interest  in  coastal  hydrooptics. 

Hydrosol  particles  in  this  model  are  perceived  as  small 
potential  holes  in  a  Brownian  motion.  They  are  characterized 
by  the  size  distribution  function:  (p{a)  (J^  (p{a)da  =  1).  Such 
a  model  allows  one  to  derive  the  following  interaction 
Hamiltonian  between  photons  and  thermal  density 


fluctuations  [6]: 


dt  ot 


1.2.3,  (1) 


here  is  a  photon  field  operator  which  corresponds  to  a 
vector  potential  of  light  wave,  r  is  a  coordinate,  t  is  time, 
index  a  denotes  vector’s  component,  c  is  the  speed  of  light, 


^  Cy  5e 


kj 


~yCy,  7  =  l0'^l0^ 


(2) 


Cy  is  the  volume  concentration  of  hydrosol,  is  the 
difference  between  average  DP  of  the  hydrosol  and  DP  of  pure 
water,  is  the  Boltzmann’s  constant,  p  is  the  water 
density,  and  u  is  the  average  velocity  of  a  Brownian 
movement  in  water.  Repeating  indices  everywhere  in  this 
article  imply  summation.  The  operator  of  thermal  density 
fluctuations  i/>  is  determined  by  the  following  equation  [6]: 


=  —  A(r  -  r')  4)(r',  t)  dr', 

^  4k  ^  ■’  (3) 

A(r)  =  l"(pia)  da  j  [4;r|~  r^dr^°(p{a)  da] , 

here  a  =  ]^(p{a)ada  is  a  mean  radius  of  hydrosol  particles,  T 

is  the  absolute  temperature. 


The  radiative  transfer  theory,  predominantly  used  in  ocean 
optics,  cannot  explain  nonlinear  dependence  of  inherent 
optical  properties  on  concentration  of  scatterers.  In  order  to 
investigate  this  problem  we  should  start  from  the  Maxwell  s 
equations  in  stochastically  scattering  medium.  The 


THE  GREEN'S  FUNCTION 
Let  us  write  a  Green's  function  of  photons  propagating  in 

a  nonscattering  medium  with  the  dielectric  permittivity  Eq.  It 

can  be  represented  as  a  sum  of  the  transverse  and  longitidinal 
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components  [6].  In  the  energetic-momentum  representation 
the  Green's  function  has  the  following  form: 


D%\o),k)  =  Z)o"-(co,k)(<5„^  -n„np+  D'(fi),k)  (4) 


where 


D^{(o,k)  =  - 


4n 


,(o^  Ic^ 


Dl,(co,k)  = 


4n:c 


(5) 


n^  is  the  component  of  the  unity  vector  in  the  direction  of 
is  the  Kroneker’s  symbol  or  the  unity  tensor. 

The  Green's  function  of  photon  field,  that  includes 
multiple  scattering  due  to  the  interaction  with  thermal  fluc¬ 
tuations  described  by  the  Hamiltonian  (1),  can  be  written  as: 
Dap{co,k)=  D" (a, k)(S^p  -n^np)+ D‘{co,k)  n^n^,  (6) 


D'"(ft),k)  =  - 


An 


co'^/c^-k^ 


D‘{(o,k)  = 


4kc 


(7) 


here  and  E^  are,  correspondingly,  the  transverse  and  the 
longitidinal  components  of  the  following  dielectric 
permittivity  tensor: 

(“>  (<«- k) (5^  £'(00, k)  .  (8) 

It  is  clear  from  Eqns  (6)-(8)  that  the  problem  of  finding 
the  dielectric  permittivity  is  equivalent  to  the  problem  of 
finding  the  Green's  function. 

Let  us  calculate  the  multiple  scattering  Green's  function 
and  the  corresponding  dielectric  permittivity.  As  a  starting 
zero  approximation  let  us  take  the  Green's  function  (4)  which 
corresponds  to  the  clear  water.  The  DP  of  the  clear  water 
depends  only  on  the  circular  frequency  co : 

£°P  (00)  =  £„  (00)  <5^  S  £„  (00)  (5„^  -  n^np )  +  fg  (oo)  .  (9) 

It  means  that  Eqn.  (9)  takes  into  account  only  temporal 
dispersion  that  is  determined  by  the  processes  of  absorption 
and  emission  of  photons  by  the  water  molecules. 

The  transverse  and  longitidinal  components  of  the 
seawater  DP  that  take  into  account  processes  of  multiple 
scattering  on  the  hydrosol  particles  can  be  expressed  as: 

E^  (CO,  k)  =  Eq  (co)  +  5  (co,  k) , 

£'(m,k)  =EQ(co)-hdEfco,k),  (10) 

=£'(cu)  +  /£''(m). 

here  E^  and  E'f  are  the  real  and  the  imaginary  parts  of  the 
dielectric  permittivity  of  pure  water. 

In  order  to  calculate  5  E^  and  S  E^  let  us  carry  out  the 
standard  procedure  [4,  6]  to  calculate  corrections  to  the  Green’s 
function  starting  from  the  Hamiltonian  given  by  Eqn.  (1). 

The  Green's  function  of  the  photons  in  turbid  water  is 
expressed  through  the  following  Dyson’s  equation: 

Dap  .  k)  =  D^ap  ,  k)  +  («„ ,  k)  (<w„ ,  k)  (ffl„ ,  k) , 

-  k)  = J  A'' (q)  (t«„ ,  k  -  q)  (fi)„ ,  k  -  q ,  k) , 

h  =  a^a)^J(2n)\  (11) 

where  is  the  polarization  operator,  and  is  the  total 

apex  part  that  corresponds  to  the  sum  of  all  orders  of  light 
scattering.  The  polarization  operator  is  linked  with  the 
dielectric  permittivity  tensor  by  the  equation: 

5  (co„ ,  k)  =  {An/ col  )n^pi(o„,k).  (12) 


By  resolving  the  Dyson’s  equation  (11)  in  respect  to  we 

have  the  following  integral  equation  for  the  components  of 
the  apex  part: 

r^(k-p,k)  =  ^5„^+;2j£(''qA^(q)r<^(k-p,k-p-q)x  (13) 

£>;,;,(k-p-q)r„^,(k-p-q,k-q)D„,,(k-q)r^(k-q,k). 

In  solving  Eqn.  (13)  let  us  restrict  ourselves  to  the  scattering 
on  hydrosols  with  the  large  average  size  of  particles, 
2  a  »  A .  In  this  case  p<q<l/d«k.By  representing  the 
apex  part  as  the  sum  of  transverse  and  longitudinal  parts 
similar  to  Eqn.  (8),  integrating  over  q,  and  analytically 
expanding  to  the  real  frequency  axis,  we  obtain  the  following 
equations  for  the  components  of  the  apex  part: 

\3  2  /^/\3 


a^co"^ 


2;rW^  ^  37r^(e0^^  >  ' 


(14) 


THE  DIELECTRIC  PERMITTIVITY 

Let  us  derive  equations  that  link  corrections  to  the 
dielectric  permittivity  with  the  apex  part.  Using  Eqns.  (9)- 
(12),  integrating  over  q,  and  analytically  transferring  to  the 
real  frequency  axis,  we  have  the  following  equations: 


E^=E^-hig 


aco\co\ 


r",  £'=£„- -4yr'.  (15) 

I6nkc  6n:  E 

In  order  to  have  a  complete  system  of  equations  let  us  add  the 
dispersion  relation  taken  from  the  Maxwell's  equations: 

(16) 

Now  we  have  the  complete  system  of  five  equations  (14)-(16) 
for  the  five  complex  parameters:  k,  T^,  ,  and  In 

the  general  case  of  arbitrary  values  of  g  (which  correspond  to 
the  arbitrary  concentrations)  this  system  of  equations  has  no 
analytical  solutions.  For  the  value  of  the  parameter 

g<A^e'j3,  Cy  <2  10-^  (or  C^<2  g/m^),  (16) 

we  have  the  approximate  solution  for  the  longitudinal 
component  of  the  dielectric  permittivity  tensor:  (17) 

’1 


E  -Ec 


1-  — sinl 


-sin 


9c  '' 


CydEj^ 


2d^pu^ 


Here  pCy  is  the  concentration  of  hydrosol  particles  in 
conventional  units  (g/m^). 

According  to  Eqn.  (17)  E^  is  determined  only  by  the 
water  properties.  As  it  is  seen  from  the  Eqns.  (14)-(16)  this  is 
also  true  for  any  water  turbidity.  For  the  calculation  of  the 
transversal  part  of  the  dielectric  permittivity  let  us  consider 
only  the  case  of  a  weak  attenuation  of  light.  Let  us  write  the 
expression  for  the  absolute  value  of  the  photon  wave  vector: 

k  =  k^  +  i(a-^  b)  /  2 ,  (18) 

here  kQ  is  the  real  part  of  the  wave  vector,  a  is  the 
absorption  coefficient,  and  b  is  the  scattering  coefficient.  The 
condition  of  weak  attenuation  a-\-b«kQ  is  equivalent  to  the 
following  conditions: 

|(5£"'| « £',  or  Cy<2- 10'^  (c^  < 2  g/m^).  (19) 

In  this  case  the  solution  of  Eqns.  (14)-(16)  for  the  transversal 
part  of  the  dielectric  permittivity  is: 
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5e’'  = 
where 


.16;r 
'  V3 


£o9c  Sign  ft)  • 


3Cy6£fl 

An 


6kJ 


apu^ 


4sin[(l/3)sin  '  cK<i. 

Ve+y¥e> 


ykoC^  j 


(20) 

(21) 


SCATTERING  COEFFICIENT 


The  light  scattering  coefficient  b  is  calculated  by  the 
substitution  of  Eqn.  (20)  into  Eqns.  (16)  and  (18).  It  has  the 
following  form: 


where 


b  = 


K  Cy5e„ 
&e'oX 


4sin[  —  sin  ‘  ^ 


,^<1, 


SCySe^  jeksT 

IEqX  ^apu^ 


(22) 

(23) 


here  X  is  the  wavelength  of  light  in  vacuum. 

If  (^  « 1  Eqns.  (22)-(23)  can  be  further  simplified: 


AE'^pd^u^?^  ’ 


Cy  <  10-=(orC^  <  10-'  g/m').  (24) 


Equation  (23)  coincides  with  the  formula  for  the  scattering 
coefficient  derived  in  Ref.  [6]  for  the  case  of  X«2a, 


In  the  general  case  of  the  arbitrary  values  of  the  parameter 
g  (or  arbitrary  concentrations  Cy  or  C^)  the  system  of  Eqns. 

(14)-(16)  has  no  solutions  in  analytic  form.  However  it  is 
easy  to  solve  that  system  numerically.  Figure  1  shows  the 
dependence  of  the  scattering  coefficient  of  seawater  as  a 
function  of  the  hydrosol  particles  concentration.  The  absolute 
values  of  b  are  not  shown  because  they  depend  on  the  average 
radius  of  scattering  particles  a  and  the  average  density  of  the 
particles  p  that  vary  from  a  region  to  a  region.  The  idea  of 
displaying  this  picture  is  to  show  the  dependence  of  the 
scattering  coefficient  b  on  concentration.  At  very  small 
concentrations  of  particles  (C^  <  0.03  mg /m^)  the  dependence 
b{Cp)  is  quadratic.  At  the  concentrations  typical  to  the  open 


0.01  0.1  1  10  100 
Concentration  Cp ,  g/m^ 


Figure  1 ,  The  scattering  coefficient  of  seawater  filled  with  the 
hydrosol  particles  as  a  function  of  the  hydrosol 
particles  concentration. 


ocean  the  dependence  is  linear.  At  higher  concentrations 
(5  <  <  50mg/m^)  the  dependence  shown  in  Fig.  1 

reminds  the  exprerimental  one  proposed  in  Ref  [1].  At  even 
larger  concentrations  (C^  >50  mg/ m^)  it  reaches  saturation. 

CONCLUSION 

It  is  shown  that  the  approach  based  on  the  Maxwell's 
theory  in  a  stochastically  scattering  medium  can  be 
productively  used  in  ocean  optics.  From  the  solutions  of  the 
Dyson’s  equation  the  dielectric  permittivity  tensor  of  seawater 
with  the  scattering  particles  is  found.  From  the  equations  for 
the  dielectric  permittivity  a  nonlinear  dependence  of  the 
scattering  coefficient  of  water  is  derived.  It  is  shown  that  the 
scattering  coefficient  of  seawater  quadratically  depends  on  the 
concentration  when  the  concentration  of  scattering  particles  is 
very  small.  It  is  linear  to  the  concentration  of  scatterers  at 
values  typical  to  the  open  ocean.  At  the  concentrations 
typical  to  the  coastal  waters  the  concentrational  dependence 
weakens  and  reaches  saturation  at  values  higher  than 
lOOg/m^. 
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